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Fault, location using measurements of current and voltage from 



one end of a line 



TECHNICAL FIELD. 
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The present invention is concerned with a fault location 
technique for a section of a power transmission line utilizing 



measurements of current and voltage made at terminals located 
at one end of the section of the power line. 



10 BACKGROUND ART 

Several methods and approaches for fault location in high 
voltage power transmission systems, and power distribution 
systems, have been developed and employed. One approach has 
been to use voltage /current transducers located at terminals 
15 located at each of two ends of a section of the power line to 
be monitored. Inductive current transformers are used to 
provide a measurement of instantaneous current in a 
transmission line. 

20 US 4,559,491 entitled Method and device for locating a fault 
point on a three-phase transmission line, describes a one-end 
fault location (FL) algorithm. High accuracy of fault location 
using a fault locator device at one end of a line is achieved 
by taking into account the actual distribution of a fault 

25 current in the transmission network. This algorithm has been 
successfully implemented into a product in 1982 and is in 
operation with single and parallel transmission lines in many 
countries around the world. However, for certain fault 
conditions it is difficult to obtain accurate pre-fault. 

3 0 quantities, such as pre-fault currents, in order to calculate 
an estimate for voltage drop across the fault path. Also, a 
disadvantage of using phase voltages and currents and zero 
sequence components of currents is that it is relatively 
difficult using these values to compensate for shunt 
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disadvantage of using phase voltages and currents and zero 
sequence components of currents is that it is relatively 
difficult using these values to compensate for shunt 
capacitance effects. In addition, the fault locator method 
described is not suitable for single and parallel line 
sections which have an extra link across the ends of the 
sections . 



SUMMARY OF THE INVENTION 

The aim of the present invention is to remedy one or more of 
the above mentioned problems . 

This is obtained by a method characterised by claim 1. 
Specific features of the present invention are characterised 
by the appending claims . 

In one aspect of the invention, a method comprising a new 
formulation of a one-end fault locator algorithm has been 
proposed. The uniform description of the transmission network 
in terms of symmetrical components as well as the generalized 
models of fault loops and faults have been applied. The 
resulting advantages include the algorithm can be used for 
locating faults in typical single and parallel transmission 
lines, and, in addition, fault location may also be carried 
out for both single and parallel lines with an extra link 
between the line ends. Another advantage is that a procedure 
for calculation of a distance to fault is in the form of a 
compact quadratic equation with the coefficients dependent on 
a fault type, acquired measurements and impedance data for the 
transmission network. Another advantage of the invention is 
that optimal estimation of the voltage drop across a fault 
path is applied, which has the result that the pre-fault 
currents in case of single phase -to -ground faults and phase- 
to-phase faults are no longer required. 
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In an embodiment, compensation for shunt capacitances is 
facilitated by means of the use of the notation of symmetrical 
components. The distributed long line model of the line has 
been applied for that. The compensation is performed 
individually for all the sequences . The currents for 
particular sequences are compensated for the shunt currents 
and then the fault loop compensated current is composed*. In 
another embodiment improved accuracy has been obtained by 
means of an option to measure the source impedance at the 
remote end instead of using a representative value. The source 
impedance measured- at the remote end may be considered as sent 
to the fault locator by using a simple communication means. 

In another embodiment, a method for one end fault location for 
parallel lines to locate single phase-to-ground faults is 
described under a plurality of conditions. In another further 
embodiment a method is described for one end fault location 
with standard availability of the measured signals for ground 
faults including both single phase- to -ground faults and phase- 
to -phase- to -ground faults. 

In another aspect of the invention, a fault locator device for 
carrying out the method of the invention is characterised by 
claim 18. Specific features of the inventive fault locator 
device are characterised by the appending claims. 

In another aspect of the invention a computer program is 
described for carrying out the method according to the 
invention. In another aspect of the invention a computer 
program product comprising a computer program for carrying out 
the method of the invention is described. In another, further 
aspect of the invention a graphical user interface is 
described for displaying a distance to a fault from one' end of 
a section of a power line. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the method and system of the 
present invention may be had by reference to the following 
5 detailed description when taken in conjunction with the 
accompanying drawings wherein: 

Figure 1 shows in a single schematic diagram a method of fault 
location in power transmission and/or distribution systems for 
10 parallel lines and single lines according to an embodiment of 
the invention; 

Figure 2a shows a schematic circuit diagram for a parallel 
transmission network for the positive sequence component in 
15 which the fault loop is marked for the cas^ of a fault locator 
installed at the terminal AA. Figures 2b, 2c show 
corresponding diagrams for the negative sequence and zero 
sequence components, respectively; 

20 Figure 3a is a schematic block diagram for obtaining and 
calculating the phasors of the symmetrical components of 
voltages and currents used for composing the fault loop 
voltage. Figure 3b shows a corresponding diagram for composing 
the fault loop current; 

25 

Figure 4 shows a circuit diagram for determining the fault 
current distribution factor for the positive sequence of a 
single line, in which diagram quantities for the negative 
sequence are shown indicated in brackets ; 

30 

Figure 5 shows a circuit diagram corresponding to Figure 4 for 
single lines for determining the fault current distribution 
factor for the positive sequence of parallel lines, in which 
quantities for the negative sequence are also shown indicated 
3 5 in brackets; 
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Figure 6 shows a schematic diagram for an embodiment of the 
invention in which source impedance measured at a remote end B 
may be communicated to a fault locator at the first end A; 

5 Figure 7 is a circuit diagram of an embodiment in which the 
shunt capacitances are taken into account, and shows a 
positive sequence circuit diagram during a first iteration; 

Figure 8 shows a negative sequence circuit diagram for taking 
10 the shunt capacitances effect into account during a first 
iteration; 

Figure 9 shows a zero sequence circuit diagram for taking the 
. shunt capacitances effect into account during a first 
15 iteration; 

Figure 10 shows a flowchart for a method for locating a fault 
in a single line according to an embodiment of the invention; 

20 Figure 11 shows a flowchart for a method for locating a fault 
in parallel lines according to an embodiment of the invention; 

Figure 12 and Figures 13a f 13b, 14, 15a and 15b show schematic 
diagrams of possible fault-types (phase-to-phase, phase- to- 
25 ground and so on) with respect to derivation of coefficients 
for Table 2 in Appendix A2 . Figure 12 shows fault types from 
a-g, and Figures 13a, 13b faults between phases a-b. Figure 14 
shows an a-b-g fault . Figures 15a and 15b show symmetrical 
faults a-b-c and a-b-c-g respectively; 

30 

Figures 16 and 17 show schematic diagrams for the derivation 
of the complex, coefficients in the fault current distribution 
factors for the positive (negative) sequence included in Table 
3. Figure 16 shows the case of a single line with an extra 
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link between the substations. Figure 17 shows the case of 
parallel lines with an extra link between the substations; 



Figure 18 shows a fault locator device and system according to 
an embodiment of the invention; 

Figure 19 shows a flowchart for a method for locating a single 
phase-to-ground fault in parallel lines in the case of • 
measurements from the healthy line being unavailable, 
according to an embodiment of the invention; 

Figure 20 shows a schematic diagram for a method of fault 
location for parallel lines with different modes of the 
healthy parallel operation; 

Figures 21a, shows a schematic equivalent circuit diagram for 
a parallel network for the incremental positive or the 
negative sequence. Figure 21b shows an equivalent circuit 
diagram for the zero sequence while both parallel lines are in 
operation. Figure 21c shows the equivalent circuit diagram for 
the zero sequence with the healthy parallel line switched off 
and grounded; 

Figure 22 shows a flowchart for a method for locating phase- 
to-phase and phase-to-ground faults in parallel lines in the 
case of providing the zero sequence currents from the healthy 
parallel line according to another embodiment of the 
invention; 

Figure 23 shows a schematic diagram fault location for 
parallel lines with standard availability of measurements 
according to another embodiment of the invention; 
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Figures 24 a,b,c show the equivalent circuit diagrams of 
parallel lines for positive, negative and zero sequence 
currents respectively. 



5 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig.l presents a schematic diagram for one-end fault location 
applied for parallel lines and for a power transmission or 
distribution system with a single line. A fault locator 1 is 
positioned at one end 2 of a single line AA-BA 3 or parallel 

10 lines AA-BA, AB-BB, 4. A fault F is shown at FA with a 

corresponding fault resistance, 5, denoted as R F . A value for 
distance to the fault d from one end 2 determined and provided 
by the fault locator 1 is indicated with the reference number 
7. Components such as parallel line AB-BB and quantities such 

15 as a parallel line value zero sequence current T^o shown with 
dotted lines are excluded when considering a single line case. 

The fault locator 1 positioned at the first end 2, or X A' end, 
is supplied with the following input signals: 

2 0 - three-phase voltages (V^) of the faulted line, 

- three-phase voltages ( of the faulted line, 

- zero sequence current (Xabo) from the healthy parallel line 
(zero sequence is not present when the single line only is 
considered) . 

25 

Fig. 2 a,b,c show circuit diagrams of a parallel transmission 
network for the positive 2a, negative 2b, and zero sequence 2c 
components. Fault loops for the sequence components 21a, 21b, 
21c are shown for the case of a fault locator installed at the 

3 0 terminal AA. An extra link 25 between the terminals A, B is 

shown. A generalized model of the fault loop considered for 
different fault types is stated as: 

^AA_p -JZilaIaajp -Rf(£fiIfi +a F ilF2 +4foIfo) = 0 (1) 
where : 
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d - unknown and sought distance to fault, 
Zlxia - positive sequence impedance of the faulted line, 
Y-AA_$> 1-aajp fault loop voltage and current composed according 
to the fault type, 
5 R F - fault resistance, 

l Fi - sequence components of the total fault current (i = 0, i 
= 1, i = 2) , 

a Ff . - weighting coefficients (TABLE 2) . 

10 Fault loop voltage and current can be expressed in terms of 
the symmetrical components of measured voltages / currents : 

V-aa_ v =*iVaai +£iY-aai +2oZaao (2) 
where : 

15 AA, AB - subscripts used for indicating measurements- 

acquired from the faulted line (AA) and from the healthy line 
(AB) , respectively, 

a Q , a x , a 2 - coefficients which are gathered in TABLE 1 (the 
Tables are arranged below at the end of the description of 
20 embodiments and derivation of these coefficients is shown in 
APPENDIX Al, also attached) . 

Zloia > £om - impedance of the faulted line and mutual coupling 
between the lines for the zero sequence, respectively, 
£om = £o " for parallel lines, 
25 a 0m =0 - for single lines. 



The phasors of symmetrical components of voltages, positive: 
Y-aai > negative: V-AA2 and zero sequence: Y-aao as well as the 
phasors of symmetrical components of currents, positive: Laa\ » 
3 0 negative: Laaz * zero sequence from the faulted line: Laao and 
zero sequence from the healthy line: l AB0 are calculated from 

8 
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the acquired measurements as shown schematically in schematic 
block diagrams Figures 3a and 3b. 



Figure 3a shows an input of instantaneous phase voltages 30a/ 
filtering stage 33a, phasors of phase voltages 31a, 
calculation of phasors of symmetrical components 33b and 
phasors of symmetrical components of voltages output at 32a. 
It may be seen from Fig 3a that acquired phase voltage 
measurements are subjected to a filter, then calculations are 
made to find the symmetrical components of the fault loop 
voltage. Figure 3b shows correspondingly stages used to find 
the symmetrical components of the fault loop current- Figure 
3b shows instantaneous phase currents and instantaneous zero 
sequence current from the healthy line 3 0b, filtering 3 3b, 
phasors of phase currents and phasor of zero sequence current 
from the healthy line 31b, calculation 34b and phasors of 
symmetrical components of currents output at 32b. 

Fault loop signals may be composed according to formulae (2)- 
(3) and TABLE 1, which is the alternative to the classic 

approach (TABLE 1A, fault loop voltage (Zaa_ p ) and current 
(Z;4A_p)' which was used in the fault locator from [1-2]. 



Voltage drop across a fault path resistance, the third term in 
(1) , can be expressed in terms of the current distribution 
factors and local measurements of currents which results in: 



ALaai , n Laa2 , _ Laao 

£F1 £F2 £F0 J 



±F2 



= 0 



(4) 



Formula (4) has been obtained from the following relations 
between the symmetrical components of a total fault current 
and measured currents : 
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t — ^Laai r _ Laxi t — Laaq 

£±F1 ~F2 —FO 

(5) 

where : 

Lfi ' Lf2' Lfo ~ symmetrical components of a total fault 
current, 

k-F\ » k.F2 ' %-fo ~ fault current distribution factors for 
particular sequence quantities, 

A/^j =Laai -Iaaiptc ' Laai ; Laaq - symmetrical components of 
currents measured in the line A at the station A (subscript 
AA) ; note that in case of the positive sequence the 
incremental quantity (post-fault current minus pre-fault 
current) is used. 

Voltage drop across the fault path, as shown in the third term 
in equation (1) , is expressed using sequence components of the 
total fault current. The weighting coefficients a F0 , a F1 , a F2 , 
can accordingly be determined by taking the boundary 
conditions for particular fault type. See TABLE 2, Alternative 
sets of the weighting coefficients for determining a voltage 
drop across the fault path resistance. Examples of derivation 
of these coefficients are contained in APPENDIX A2 . 



There is some freedom for setting the weighting coefficients. 
It is proposed to utilize this freedom firstly for avoiding 
zero sequence quantities, since the zero sequence impedance of 
a line may be considered as an unreliable parameter. This can 
be accomplished by setting a F0 =0 as shown in TABLE 2. 

Secondly, the freedom in establishing the weighting 
coefficients can be utilized for determining the preference 
for using particular quantities. The negative sequence (TABLE 
2, set I) or the positive sequence (TABLE 2, set II) can be 
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preferred as well as possibly both types of the quantities 
(TABLE 2, set III) can be used for determining the voltage 
drop across the fault path. 

The set I is recommended for further use, thus avoiding, the 
positive sequence, and thus avoiding the pre-fault positive 
sequence current, for the largest number of faults. Avoiding 
the pre-fault positive sequence current is highly desirable 
since sometimes the pre-fault currents - due to certain 
reasons - can not be recorded or registered, but may be 
contaminated by one or more the symptoms of the occurring 
fault. Moreover, the accuracy of recording the pre-fault 
currents, which are basically lower than the post-fault 
currents, is not very great. This is so since the A/D 
converters operate with less accuracy in the low range. 

Fault current distribution factors depend on the configuration 
of the transmission network, Figures 4, 5, and impedance 
parameters. Basically, all impedances for the positive and for 
the negative sequence are equal to each other and thus one 
obtains : 

k -1c -E^ah, 

~ F1 ~~ F2 ~m x — (6) 

Coefficients in a fault current distribution factor (6) for a 
single (Figure 4) and for parallel lines (Figure 5) are 
gathered in TABLE 3, Coefficients for determining a fault 
current distribution factor, (note that derivation of the 
coefficients is shown in APPENDIX A3) . 

Figure 4 shows a circuit diagram of - a single line for 
determining the fault current distribution factor 
for the positive sequence currents and with the negative 
sequence currents such as shown in brackets. Similarly 
Figure 5 shows a circuit diagram of parallel lines for 
determining the fault current distribution factor with 
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positive sequence currents wherein the negative sequence 
currents are shown in brackets . 



In Fig. 4 the extra link 45 between the terminals A, B having 
impedance for the positive sequence equal to: Z lAB can be 
considered as existing (Z {AB #c ) or as not present (Z 1AJ? -»oc) . 
In Fig. 5 the extra link 55 between the terminals A, B having 
impedance for the positive sequence equal to: Z 1AB can be 

Z Z 

considered as existing ( Z lLB8cAB = S~ lLB ~™ B ) or as not present 

Zlilb + Z-\ab 

( Z-\lb8lAb " Z lLB ) . 

Substituting (6) into (4) and adjusting a FO =0 (as in TABLE 2) 
results in: 



After multiplying both sides of (7) by: =J — and some 

rearrangements, the quadratic formula with two unknowns, d- 
[p.u.] sought fault distance from A, R F - fault resistance, is 
obtained: 




~ dZ lLA LAAjp 




(7) 



KiZ 1L d 2 + (L X Z 1L -K+Z^d -UZaaj* + RfKi 



(&fi a Iaa\ + QlfiLaat) = q 



(8) 



where : 




calculated fault loop impedance. 



Writing formula (8) in more compact form results in: 
A 2 d 2 + A x d + A 0 +A 00 R F = 0 



(8a) 



where : 



Al = 4z_Re + J A 2_lm = K\Z.1LA 

Al = A_Re + Ml_Im = L\Z\LA ~~ KlZj^jp 

Ao = A)_Re + jA)_Im = -£lZ/iA_ P 
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A .i A _ KMfi a Laai +*fiIaa2) 

A)0_Re JA)0_Im — T 

V AA 

Z AA n - ~ ~ p - calculated fault loop impedance 
— aa_p j 

£1 , £1/ Mi - coefficients gathered in TABLE 3. 

Formula (8a) can be written separately for real and imaginary 
parts : 

A 2 Ke d 2 + A l Kc d + ,V Re + 4»_Re*F =0 ( 8b) 

A 2 ^d 2 + Ai_ ta d + Vim + 4»_Im*F =0 ( 8C) 

Combining (8b) and (8c) in such the way that fault resistance 
is eliminated [that is, equation (8b) is multiplied by ,A)0_im 
and equation (8c) by i4 00 _ Re and then subtracting them] yields 
the quadratic formula for a sought fault distance: 
B 2 d 2 +S^+B 0 =0 O) 
where : 

B 2 = ^2_ReA)0_Im ~" ^2_ImA)0_Re 
B\ = A_ReA)0_Im "" A_ImA)0_Re 
#0 ~ A)_ReA)0_Im ~ A)_ImA)0_Re 

Equation (9) has two roots (d lf d 2 ) for the distance to fault: 



_ -B 1 -tIb*-4B 2 B 0 

d ^ W 2 



J -B x + ,jB?- 4B 2 B 0 
d 2 — — _ 



2B 2 



(10) 



The root which fulfils the condition (0<d<l) is selected as 
the solution for the distance to fault. 
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In another embodiment of the invention, the method of fault 
location is carried out by using a measurement of source 
impedance at the second end remote from the fault locator 1, 
instead of a representative value for the source impedance at 
5 the remote end, and by communicating that measured value to 
the local end using a communication means. Coefficients- from 
(9) are determined with the local measurements and the 
impedance data for the transmission line, the extra link 
between the line terminals and the equivalent systems at the 
10 line terminals. Impedance of the equivalent system at the 
local substation {Z lsA ) can be traced on-line with the local 
measurements. In contrast, the remote system impedance (Z isB ) 
is not measurable locally from A. Thus, the "representative" 
value of this impedance may be provided for the algorithm [1- 



The alternative solution for the single line case is shown in 
Figure 6, which shows a fault locator 1 at the first end 2 
near to system A, and another device 10 located at the remote 
20 end 3 close to system B, indicated as RD. A communication 

signal 9 is shown being sent from the 10 at the remote end to 
the fault locator 1 at the local end. 

The remote source impedance (Z lsB ) is measured by the remote 
25 device RD, 10, which may be another fault locator or any 
suitable device such as a digital relay or digital fault 
recorder, in the remote substation and the measurement 9 is 
sent via a communication channel 60. Synchronization of 
measurements at the line terminals is not required. The source 
3 0 impedance is calculated from the known relation between the 
incremental positive voltage ( AV Bl ) and the incremental 
positive sequence current ( Al Bl ) [3-4]: 



15 



2] . 




(11) 
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Similarly, the fault locator 1 calculates the local source 
impedance 

Z lsA ^-—~ (Ha) 
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In another and preferred embodiment of the invention, 
compensation is carried out for the shunt capacitance of the 
line. Compensation for shunt capacitances effect can be 
accomplished by taking into account the lumped model (only the 
longitudinal R-X L parameters are taken into account) or the 
distributed long transmission line model. The distributed long 
line model [5] as providing higher accuracy of fault location, 
has been considered here. 

The compensation for the single line is presented further. 
15 This means that when composing fault loop current (3) the term 
reflecting the mutual coupling effect disappears (a 0m =0). 
Moreover, the single subscript {A instead of AA) is used. 



Fault location procedure with compensating for shunt 
20 capacitances of a transmission line requires the following 
additional input data: 

C 1L - shunt capacitance of a whole line for the positive and 
the negative sequences (parameters of a line for the positive 
and the negative sequences are identical and thus: C 2L = C lL ) 
25 C 0L - shunt capacitance of a whole line for the zero sequence, 
I - total line length (km) , used for expressing 
impedances /capacitances of the line per km length. 

The compensation of shunt capacitances may be introduced while 
3 0 determining the voltage drop across the faulted line segment - 
the second term in the generalized fault loop model (1) . This 
requires compensating the components of the calculated 
currents for particular sequences. Thus, the original measured 
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currents : l A1 , J A2 , l A0 have to be replaced by the currents 

after the introduced compensation: LAi^comp > Lxi^comp * l-AQ^comp • 
At the same time the original fault loop voltage, the first 
term in the model (1) , is taken for a distance to fault 
calculation. As concerns determining the voltage drop across 
the fault resistance, the third term in (1), it is assumed 
here, which is a standard practice, that the effect of line 
capacitances at the location of the fault (point F) , may be 
neglected. This is justified as the impedance of the 
capacitive branch at that location is much greater than the 
fault resistance. This means that the voltage drop across the 
fault resistance is determined without taking into account the 
shunt capacitances . 

Calculating a distance to fault the following impedances 
(defined below) are taken as: 

Zl\l 8 ~ positive sequence impedance of a line with taking into 
account the distributed long line model, 
Zlql 8 " as above, but for the zero sequence. 

The compensation procedure requires iterative calculations, 
performed until the convergence is achieved (i.e. iterated 
until the location estimate ceases to differ from the previous 
estimate) . However, the studies conducted revealed that 
results of acceptable accuracy may be obtained using a fixed 
number of iterations, for example, 2-3 iterations. The 
calculated distance to a fault from a particular (say, present 
iteration) is utilized for determining the shunt current in 
the next iteration. The determined shunt current is then 
deduced from the measured current. A distance to fault, 
calculated without considering the shunt effect (10) , is taken 
as the starting value for the first iteration. 
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A way of conducting the first iteration of the compensation is 
shown in Figures 7, 8, 9 for the positive sequence, negative 
sequence and zero sequence respectively with taking into 
account the shunt capacitances effect. 

As a result of performing the first iteration for the positive 
sequence (Figure 7) the compensated current LA\_comp_\ *- s 
calculated and the last index in the subscript denotes the 
first iteration. The calculation is based on deducing the 
shunt current from the calculated positive sequence current 7 A1 
calculated from the measured phase currents - Fig. 2: 

L A \_comp_\ = Lai ~ O^v^lA^V^ ( 12 ) 

where : 

d v - distance to fault calculated without taking into account 
the shunt capacitance effect (10) , 
I - total line length (km) 



t 

^tanhl = 



anh^0.5Z; L 5i L rf v /J 
*Jo.5Z 1L B' lL d v l 



B[ L =i— 1^- - positive sequence admittance (capacitive) of a 
line per km length (S/km) 

Z\ L = - positive sequence impedance of a line per km length 

(Q/km) 

Positive sequence impedance of a faulted line segment, 
between points A and F , without taking into account the shunt 
capacitances effect and considering the lumped line model 
equals : 

d v lZ 1L (13) 

while for the considered here distributed long line model: 
dJZ^A^ (14) 
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where : 



^Z lL B lL d v l 



Thus, the positive sequence impedance of a line with taking 
into account the distributed long line model {Z 1 ^ 8 ) equals: 
Z% 8 =A sMll Z lL (15) 

As a result of performing the first iteration for the negative 
sequence (Figure 8) the compensated current L A 2_comp_\ is ' 
calculated and the last index in the subscript denotes the 
first iteration. This is based on deducing the shunt current 
from the calculated negative sequence current 7 A2 , calculated 
from the measured phase currents - Fig. 2) : 

L A2 _comp_l =lA2 ^d v lB 2L A^ 2 V A2 < 16 > 

where, taking into account that the line parameters for the 
positive and for the negative sequences are identical 
( C 2L = C lL , Z 2L — Z lL ) : 

4tanh2 = ^tanhl 

As a result of performing the first iteration for the zero 
sequence (Figure 9) the compensated current l A o_com P _\ is 
calculated, last index in the subscript denotes the first 
iteration. This calculation is based on deducing the shunt 
current from the calculated zero sequence current 7 A0 , 
calculated from the measured phase currents - Fig. 2: 

LAO _ compel = LaO ~ O-^v^OL^tanhoZAO ( 17 ) 

where : 
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AtznhO 



tanh^Jo.5Z OL B QL d v l^ 
^0.5Z 0L B 0L d v l 



S.hL = — " zero sequence admittance (capacitive) of a line 
per km length (S/kiti) 

Z' 0L ==^ - zero sequence impedance of a line per km length 
I 

(Q/km) 



Zero sequence impedance of a faulted line segment, that is 
between points A and F, without taking into account the shunt 
capacitances effect and considering the lumped line model 
10 equals: 

d v lZ 0JL (18) 

while for the considered here distributed long line model: 
dJZ 0L A shih0 (19) 



15 where: 



smh^Z 0L B 0L d v l^ 



^sinhO 



Thus, the zero sequence impedance of a line with taking into 
account the distributed long line model (Zjj^M equals: 
20 Z$ s =A sinh0 Z 0L (2 0) 

The quadratic complex formula (8) with two unknowns 

( d comp t [p.u. ] - sought fault distance, R F - fault resistance) 

after introducing the compensation (first iteration) takes the 
25 following form: 

(21) 

19 
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where : 

1 V A 

Z c °"*p- = "omplli ~ fa ult loop impedance calculated with: 
Ia_ p 

V Ap - original fault loop voltage (2), 

-A_p —S-\LA\_comp_\ + 2.2±-A2_comp_i + Q.o4-A0_comp _i ~ fault loop current 
composed of the positive (12), negative (16) and zero (17) 
sequence currents obtained after deducing the respective 
capacitive currents from the original currents. 

Writing (21) in a more compact form results in: 

^^d^0 2 +AT'- 1 d a ^_ l +Ar f - l +A^^^0 (21a) 
where : 

jcamp_ 

1 V a 

Z c ™p- = ~ - v x - fault loop impedance calculated with: 
Ia_ p 

¥-a_ p ~ original (uncompensated) fault loop voltage (2), 

jcomp_l _ j j 

La_p -£i4-Ai_comp_i + £2L A 2_comp_i+a.olAO_comp_i ~ fault loop current 

composed of the positive (12), negative (16) and zero (17) 

sequence currents obtained after deducing the respective 

capacitive currents from the original currents, 

Ki / L\, Mi - coefficients gathered in TABLE 3 . 

Formula (21a) can be written down separately for real and 

imaginary parts : 
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^ c J^ compl f + Affr^j + + AfiPfc 1 *, = 0 < 21b > 
^TSr 1 W-v j) a + K3t 1 *o^jl + ^r 1 + = 0 ( 21c) 

5 Combining (21b) and (21c) in such the way that fault 
resistance is eliminated, that is, equation (21b) is 
multiplied by A$% 1 and equation (21c) by ^fc 1 and then 
subtracting them, yields the quadratic formula for a sought 
fault distance: 

10 Br P -Hd comp _^ 2 +B^-% omp _ 1 +Br P - 1 =0 (22 ) 

where : 

Ttcomp^.1 _ a compel a compel _ AComp_l a compel 
&2 ™2_Re ^O.Im A 2_Im ^O^Re 

j^comp _\ _ pcomp _ 1 ^i^™^ * — A^°^~ * i ^Oo' 7 ^R ^ 

n compel „ ,ic<wz/?_l /* compel AComp_lAComp_l 
B 0 -^O^Re Ax)_Im ^0_lm ^O.Re 
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Equation (22) has two roots [ (^ comp ^i)i / (^m P _i)2 3 for the 
distance to fault: 



^ camp _i ^ (5 comp - 1 ) 2 — AB^ mp ~~ 1 BQ° mp - 1 

2 * 2 (23) 

(.d com p_\)2 "~ t% compel 



20 The root, which corresponds to the selected previously root 
(10) for d (uncompensated) is taken as the valid result. The 
compensation procedure requires iterative calculations, 
performed until the convergence is achieved (i.e. until the 
location estimates cease to change from the previous 

25 estimates) or as with a fixed number of iterations such as 2-3 
iterations. The calculated distance to a fault from a 
particular (say, present iteration) is utilized for 
determining the shunt current in the next iteration. 
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The method of the invention is illustrated in two flow-charts 
of the FL algorithm, Figure 10, single line and in Figure 11, 
parallel lines. 

As shown in the flowchart in Figure 10 the following 
measurements are utilized: 

- voltages from the side A from particular phases a, b, c: 
v A_a > v A b , v A c 

- currents from the side A from particular phases a, b, c: i A 
' *A_b * i>A_c 

The input data utilized at step 101, Input data and 
measurements, are as follows: 

- impedances of the line for the positive (Z lL ) and zero (Z 0L ) 
sequences , 

- impedance of the extra link 25, 45, 55 between the 
substations A, B for the positive (negative) sequences (Z 1AB ) 

- source impedances for the positive (negative) sequences 
(Z.isA> Zlub) : the representative values or the measured values 
are used, and a communication means is used for sending the 
measured remote source impedance as previously described, 

- information on the fault type (from the protective relay) . 

The measured fault quantities (voltages and currents) undergo 
adaptive filtering at step 104, Adaptive filtering of phase 
quantities, aimed at rejecting the dc components from currents 
and the transients induced by Capacitive Voltage Transformers 
(CVTs) . 

In the next step the symmetrical components of voltages and 
currents are calculated, step 105, which is equivalent to that 
as shown in Figures 3a, 3b. The fault loop signals are 
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composed: the fault loop voltage as in (2) , while the fault 

loop current as in (3) - but with taking a 0m =0. 



The distance to fault without taking into account the shunt 
capacitances effect (d) is calculated at step 106 by solving 
the quadratic formula (9) . The solution of (9) is presented in 
(10) . 

The result obtained without taking into account the shunt 
capacitances effect d following 106 is treated as the starting 
value for performing the compensation for shunt capacitances . 
The distributed long line model is applied for the 
compensation . 

The following additional data is required for the calculating 
the compensation for shunt capacitance, step 107: 

- positive sequence capacitance of the line ( C 1L ) 

- zero sequence capacitance of the line (C 0L ) 

- line length (2) , which is used to express line impedances / 
capacitances per km length. 

The first iteration of the compensation leads to the quadratic 
equation (22), which is solved in (23). The next iterations 
are performed analogously. Iterative calculations are 
performed until the convergence is achieved or a fixed number 
of iterations, i.e. 2-3 iterations, may be made. The 
calculated distance to a fault from a particular (say, present 
iteration) is utilized for determining the shunt current in 
the next iteration. After completing the iterative 
calculations the distance to fault d comp is obtained. 

As shown in the flowchart in Figure 11 for parallel lines the 
following measurements are utilized: 

- voltages from the side A and line LA from particular phases 

a, b, C: v AA_a ' V AA_b ' y.AA__c 
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- currents from the side A and line LA from particular phases 
a, b, c: i AA _ a , ^ ^ , i M c 

- zero sequence current from the healthy parallel line LB: i AB0 

5 The input data utilized in step 111, Input data and 
measurements, are as follows: 

- impedances of the faulted line for the positive (Z 1M ) and 
zero ( Zlqla ) sequences 

- impedance of the healthy line for the positive (negative) 
10 s equenc e ( Z iLB ) 

- impedance of the extra link between the substations A, B for 
the positive (negative) sequences 

(Z 1AB ) 

- zero sequence impedance for the mutual coupling (Z 0m ) 
15 - representative values of the source impedances for the 

positive (negative) sequences 

( Zl\sA ' Zl\sB ) 

- information on the fault type is obtained from the 
protective relay. 

20 The measured fault quantities, the voltages and currents, 
undergo adaptive filtering at step 114 for the purpose of 
rejecting the dc components from currents and the transients 
induced by Capacitive Voltage Transformers (CVTs) . 

In the next step 115 the symmetrical components of voltages 
25 and currents are calculated as shown in Figures 3a, 3b. The 

fault loop signals are composed: the fault loop voltage as in 
(2), while the fault loop current as in (3) - but with taking 

The distance to fault without taking into account the shunt 
3 0 capacitances effect (d) is calculated in step 116 by solving 
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the quadratic formula (9) . The solution of (9) is presented in 
(10) . 

The result obtained without taking into account the shunt 
capacitances effect (d) is treated as the starting value for 
5 performing the compensation for shunt capacitances . The 

distributed long line model is applied for the compensation. 

The following additional data for the faulted line is required 
for the compensation for shunt capacitances step 117: 

- positive sequence capacitance of the line ( C lL ) 
10 - zero sequence capacitance of the line ( C 0L ) 

- line length (1) , which is used to express line impedances / 
capacitances per km length. 

In the case of a single line, the compensation is performed 
15 analogously. Iterative calculations are performed until the 
convergence is achieved or by using a fixed number of 
iterations, e.g. 2-3 iterations. The calculated distance to a 
fault from a particular iteration, for example the present 
iteration, is utilized for determining the shunt current in 
20 the next iteration. After completing the iterative 

calculations the distance to fault d comp is obtained. 

Figure 18 shows an embodiment of a device for determining the 
distance from one end, here shown as end A of a section of 
power transmission or distribution line A-B, to a fault F on 

25 the power line according to the method of the invention 

described. The fault locator device 1 receives measurements 
from measuring devices located at one end A such as current 
measuring means 14, and voltage measurements from voltage 
measurement means 11. The fault locator device may comprise 

30 measurement value converters, members for treatment of the 

calculating algorithms of the method, indicating means for the 
calculated distance to fault and a printer or connection to a 
printer or facsimile machine or similar for printout of the 
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calculated fault. In a preferred embodiment of the device the 
fault locator comprises computer program means for providing a 
display of the information provided by the method of the 
invention, such as distance to a fault d or d comp on a terminal 
5 which may be remote from the location of the line and/or the 
fault locator. Preferably the computer program means receives 
information from the fault locator and makes it available to 
provide information for a display of a computer such that an 
operator or engineer may see a value for the calculated 
10 distance to a fault displayed. The value may be displayed 
relative to a schematic display of the line or network in 
which the fault has taken place. 

In the embodiment shown, measuring device 14 for continuous 

15 measurement of all the phase currents, and measuring device 11 
for measurement of voltages are arranged at one end, station 
A. Optionally, measuring devices such as 15, 13 may also be 
arranged at station B but they are not necessary to practise 
the invention. The measured values such as: the three-phase 

2 0 voltages (Vaa) of the faulted line, three-phase voltages (Iaa) 
of the faulted line and zero sequence current {Iabq) from the 
healthy parallel line (note that zero sequence is not present 
when the single line only is considered) , and a value 
representative of the source impedance at B, Z lsB as a are all 

25 passed to a calculating unit comprised in fault locator 1, 

filtered such as described in relation to figure 3a, 3b, and 
stored in memory means . The calculating unit is provided with 
the calculating algorithms described, and programmed for the 
processes needed for calculating the distance to fault. 

30 Optionally, the source impedance for the remote end, Z lsB may 
be measured by the remote device RD # 10, and the information 
sent via a high speed communication means 60 to the fault 
locator at A. In some applications it will be preferable to 
use a measured value sent from' B instead of a representative 

35 value stored at A. It may be seen in Figure 18 that current 
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measuring means 15 and voltage measuring means 13 at remote 
end B may provide a RD 10, a fault locator, or any suitable 
device with measurements to calculate the remote source 
impedance 

5 

The calculating unit of fault locator 1 is provided with pre- 
fault phase currents and also with known values such as shunt 
capacitances and the impedances of the line. In respect of the 
occurrence of a fault, information regarding the type of 

10 fault, phase-to-phase, phase-to-ground etc., may be supplied 
to the calculating unit of the fault locator. When the 
calculating unit has determined the distance to fault, it is 
displayed on the device and/or sent to display means which may 
be remotely located. A printout or fax of the result may also 

15 be provided. In addition to signalling the fault distance, the 
device can produce reports, in which are recorded measured 
values of the currents of both lines, voltages, type of fault 
and other measured and/or calculated information associated 
with a given fault at a distance. 

20 

The method and a fault locator device according to any 
embodiment of the invention may be used to determine distance 
to a fault on a section of power transmission line. The 
present invention may also be used to determine a distance to 
25 a fault on a section of a power distribution line, or any 
other line or bus arranged for any of generation, 
transmission, distribution, control or consumption of 
electrical power. 



3 0 The fault locator device and system may comprise filters for 

filtering the signals, converters for sampling the signals and 
one or more micro computers. The microprocessor (or 
processors) comprises a central processing unit CPU performing 
the steps of the method according to the invention. This is 

3 5 performed with the aid of a computer program, which is stored 
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in the program memory- It is to be understood that the 
computer program may also be run on one or more general 
purpose industrial computers or microprocessors instead of a 
specially adapted computer. 

The computer program comprises computer program code elements 
or software code portions that make the computer perform the 
method using equations, algorithms, data and calculations 
previously described. A part of the program may be stored in a 
processor as above, but also in a ROM, RAM, PROM or EPROM chip 
or similar. The program in part or in whole may also be stored 
on, or in, other suitable computer readable medium such as a 
magnetic disk, CD-ROM or DVD disk, hard disk, magneto-optical 
memory storage means, in volatile memory, in flash memory, as 
firmware, or stored on a data server. 

A computer program according to the invention may be stored at 
least in part on different mediums that are computer readable. 
Archive copies may be stored on standard magnetic disks-, hard 
drives, CD or DVD disks, or magnetic tape.* The databases and 
libraries are stored preferably on one or more local or remote 
data servers, but the computer program and/or computer program 
product may, for example at different times, be stored in any 
of; a volatile Random Access memory (RAM) of a computer or 
processor, a hard drive, an optical or magneto-optical drive, 
or in a type of non-volatile memory such as a ROM, PROM, or 
EPROM device. The computer program may also be arranged in 
part as a distributed application capable of running on 
several different computers or computer systems at more or 
less the same time. 

In another preferred embodiment, the fault locator may be used 
with parallel lines to locate single phase-to-ground faults 
(a-g, b-g, c-g faults) in case of unavailability of 
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measurements from the healthy parallel line. Two modes of the 
healthy line operation are taken into account: 

- the healthy line being in operation, 

- the healthy line switched-of f and grounded at both the ends . 

5 Figure 19 shows the flow-chart of the algorithm for locating 
faults in parallel lines under unavailability of measurements 
from the healthy parallel line. The unavailable, here healthy 
line zero sequence current, is required for reflecting the 
mutual coupling effect under single phase-to-ground faults (a- 
10 g, b-g, c-g faults) . The unavailable current is thus 

estimated. The other faults can be located with the standard 
fault location algorithm (such as the algorithm from reference 
[1]) . 

The sequence of computations for the presented one-end fault 
15 location algorithm is as follows. 

As shown in the flowchart in Figure 19 the following 
measurement s are utilized: 

- voltages from the side A and line LA (faulted) from 
particular phases : v >u _ a , b , v AA _ c 

2 0 - currents from the side A and line LA (faulted) from 

particular phases : i M a , < 1 aa_ c 

The utilized input data are as follows: 

- impedances of the faulted line for the positive (Z 1£A ) and 
zero ( Z 0LA ) sequences 

25 - impedance of the healthy line for the zero sequence (Z QLB ) 

- zero sequence impedance for the mutual coupling (Z 0m ) 

- information on the fault type (from the protective relay) 

- mode of the healthy line operation: in operation or 
switched-of f and grounded at both ends) . 

3 0 The measured fault quantities (voltages and currents) undergo 

adaptive filtering aimed at rejecting the dc components from 
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currents and the transients induced by Capacitive Voltage 

Transformers (CVTs) . 

The following equations are used in the method of the present 
parallel line embodiment. In addition to the algorithm (1) 
described above .for estimating distance to fault (d [pu] ) by 
considering the Kirchhoff's voltage law for the fault loop as 
seen from the locator installation point: 

( a^AI d 



¥-AA_p -^ZliLaLaa^ " R F 



<±.F1^±AAI j <±F2L 
k-Fl 



QlfiLaki ! QfqLaao |_q (24) 

£~F2 ^F0 ) 



the fault loop voltage ( Y-AAjp ) anc ^ current ( LaAjp ) can be 
10 expressed in terms of the symmetrical quantities: 

¥-AA_p ="l¥.AAl +aiY-AA2 + ZqVaAO (2 ) 

and 



, SL Z 

Laa_p = Laa_p + ~ m I abo (25) 
where : 



15 /£ p ^a.lAAX +<k2lAAi +a 0 ^ L L A A0 (25a) 

—\LA 

is the fault loop current without compensating for the mutual 
coupling effect (i.e. composed as for the single line - 
superscript SL) , 

£ i ' ^2' £o"~ complex coefficients gathered in TABLE 1 (the 
20 derivation as in APPENDIX Al) , 

Y-aai ' Y-AA2* ¥-AA0 ~ positive, negative and zero sequence of 
measured voltages, 

Lam ' Laa2> Laaq " positive, negative and zero sequence current 
from the faulted line LA, 

25 Labo " unavailable zero sequence current from the healthy 
parallel line LB (to be estimated) , 
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Z 1M , Zlqla " positive and zero sequence impedance the whole 

line LA, 

Z Qm - zero sequence impedance for mutual coupling between the 
lines LA and LB, 

5 R F - unknown fault resistance. 

In the next step the symmetrical components of voltages and 
currents are calculated as shown in Figures 3a, 3b. The fault 
loop signals are composed: the fault loop voltage as in (2), 
while the fault loop current as in (25). The formulae (2) -(25) 
10 present the fault loop signals expressed in terms of the 

symmetrical components of the measured signals. However, one 
can use the classic way for composing the fault loop signals. 



The presented method covers single phase- to ground faults (a- 
15 g, b-g, c-g faults) . The other remaining faults have to be 

located with the fault location algorithms described above or 
standard fault location algorithm [1] . Distance to fault (d) 
for the considered here single phase-to ground faults is 
calculated by solving the quadratic formula for a sought 
20 distance to fault (26) . Equation (26) is the same as equation 
(10) except that the values for B X ,B 29 B 3 are different to the 
values determined in (10). The solution gives two roots: 



2B. 



2 



-b x + 4bI-ab z b 0 
d ^ w 2 



(26) 



(as above, the root which fulfils the condition (0<d<l) is 
25 selected as the solution for the distance to fault) . One has 
to substitute the following into (26) : 
B 2 = real(A 2 ) imagCAoo ) - real(A 00 ) imag(A 2 ) 
B x = realCd! ) imag(A 00 ) - real(A 00 ) imagCAj ) 
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B Q = real(A 0 ) imag(A 00 ) - real(A 00 ) imag(A 0 ) 

where : 

—2 -"^lM^llM.p Z AllAAO 

1-0 i-0 

SL Z 

400 =-M.l«F2iAA2 

The recommended SET of the coefficients & F1 , b F2 are taken from 
TABLE 4 and the recommended SET of the coefficients a Fl , a F2 , 
a WCi from TABLE 5 . 



Fault loop voltage in this embodiment is found from the TABLE 
below 



Fault loop voltage composed 
in terms of symmetrical 
components 


Fault loop voltage composed 
as in the classic approach 


a-g fault: a 2 = « 2 = a 0 = 1 

2 

b-g fault: «j = a , a 2 = a , a 0 = 1 
c-gr fault: #i = a , a 2 =a , a 0 =l 
a = exp(j27T/3) , j = V^T 


a-g fault: V M-P =E^ a 
a-gr fault: E^ -P =Zm^ 
a-gr fault: V P = ^„c 



Fault loop current 7 i4A _ p composed as for the single line is 
found from the TABLE below 
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SL 

Iaa_p i n terms of symmetrical 


p as i n the classic approach 


components 




jSL, j j ^ OLA r 
±AA_p &\LaAI +&2LAA2 + *0 ~ i-AAO 




a-gr fault: a x = a 2 = a 0 = 1 


a -gr fault: Laa.p = Laa.c + ^oLaao 


£>-gr fault : a x = a 2 , a 2 =a , a 0 = 1 


SL. 

a-gr fault: Z A4 _ p = / /U _ a + £ 0 Lwo 


c-gr fault: a j = a , a 2 ~g} , a 0 =l 


Of 

a-gr fault: Iaa_ p = L M _a + tolAAo 


a = exp(j2^73) , j = V^T 


where: fc 0 = ^0M~^iM 





The complex coefficients dependent on the mode of the healthy 
5 parallel line operation: 

a) healthy line LB is in operation: 

—OLA ~0/» 

o Mi = z 1M z 1LB +z 11A <z lsA +z lsB )+z 1LB (z lsA +z lsB ) 



b) healthy line I/B is switched-off and grounded: 




^1 = £lLA 

M!=Z lsA +Z lsA + Z ilA 
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In another embodiment, which is presented here, the method is 

applied for the standard availability of the measured signals 

and is only valid for ground faults, including both: 

- single phase- to-ground faults 

- phase-to-phase-to-ground faults. 

Fault location procedures obtained for these faults under the 
standard availability of the fault locator input signals are 
extremely simple and compact. Distance to a fault is 
calculated with a first order formula. 

Figure 22 shows a flow-chart of the developed algorithm for 
locating ground faults in parallel transmission lines. 
The sequence of computations for the presented one-end fault 
locator is as follows. As shown in the flowchart in Figure 
22. B3 the following measurements are utilized: 

- voltages from the side A and line LA from particular phases 
a, h, c: v AA _ fl , VAA^b > ^AA_c 

- currents from the side A and line LA from particular phases 
a, b, C: lAA^a ' l AA_b ' *AA_c 

- zero sequence current from the healthy parallel line LBi 
The utilized input data are as follows: 

- impedances of the faulted line for the positive ( Z lLA ) and 
zero (Z 0LA ) sequences, 

- impedance of the healthy line for the zero sequence (Z QLB ) 

- zero sequence impedance for the mutual coupling (Z 0m ) 

~ information on the fault type from the protective relay. 

The measured fault quantities (voltages and currents) undergo 
adaptive filtering aimed at rejecting the dc components from 
currents and the transients induced by Capacitive Voltage 
Transformers (CVTs) , preferably as described and shown in 
relation to Figures 3a, 3b. 
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A generalized model of the fault loop, used for deriving the 

algorithm of the present embodiment, is stated as follows: 

¥-AA_p ~ d Z lLA LAA_p ~ R F (&F\Lf\ + "fiIf! + ^FoIfO ) = 0 < 1 > 

where : 

5 d - unknown and sought distance to fault, 

Z lLA - positive sequence impedance of the faulted line, 
V_aa p> Laa_p ~ fault loop voltage and current composed according 
to the fault type, 
R F - fault resistance, 
10 Lfi " sequence components of the total fault current (i=0 - 
zero sequence, 1=1 positive sequence, 

i=2 - negative sequence) , 
a Fi - weighting coefficients (TABLE 2) . 

Fault loop voltage and current can be expressed as in classic 
15 distance protection technique or, as in this document, in 
terms of the local measurements and with using the 
coefficients {a 0 , a x , a 2 ) which are gathered in TABLE 1 
(derivation of the coefficients is contained in APPENDIX 
APP . 1 ) : 

2 0 V Mp = AA1 + a 2 VAA 2 + a Q VAAo < 2 ) 

LAA p = dllAAl + <LllAA2 + <*0 f^-LuO + &0 ^Lab* < 3 ) 

where : 

AA, AB - subscripts used for indicating measurements acquired 
from the faulted line (AA) and from the healthy line (AB) , 
25 respectively, 

Z olA , Z 0m - impedance of the faulted line and mutual coupling 
between the lines for the zero sequence, respectively. 

In the next step the symmetrical components of voltages and 
currents are calculated as shown in Figure 3a, 3b. The fault 
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loop signals are composed: the fault loop voltage as in (2), 
while the fault loop current as in (3). The formulae (2) -(3) 
present the fault loop signals expressed in terms of the 
symmetrical components of the measured signals. However/ one 
5 may instead use the classic way for composing the fault loop 
signals, as shown in APPENDIX Al . 

The presented method covers single phase- to ground faults (a- 
g, b-g, c-g faults) and phase- to -phase- to -ground faults (a-b- 
g, b-c-g, c-a-g faults), thus, the faults for which the 
10 highest fault resistance can be expected- The other remaining 
faults have to be located with the fault location algorithms 
described above or a standard fault location algorithm, such 
as for example the fault locator from reference [1]). 



15 



Distance to fault (d) for single phase- to -ground faults is 
calculated as follows : 



a-g fault: 

imag{{Z llA lAA_ v )Y$aAAo ' EoLaboW) 
where : 



(27a) 



Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as in the classic approach 


VAA_p =£lZAAl +a2¥-AA2 +ao¥-AAO 




Laa_p ~&iLaai +&iIaa2 + 

i n ~0LA r , _ Zom r 

+ ry LaAO + £0 ry I ABO 

=-\LA £L\LA 


^A4_p ~^AA_a +koi-AAO + ^0m^AB0 


where : a x = a 2 = Q 0 = 1 


where : k 0 - ~ 0M ~~ lLA , k 0m ~ 



Z —Z 

Eo = ~ 0LB ~Z° m (for symmetrical lines: P 0 =1 ) 



20 



b-g fault: 
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d _ imagiV^ ^aAAO -£oIajh>)]*> 


(27b) 




where : 




Fault loop signals composed 


Fault loop signals composed 


in terms of symmetrical 


as in the classic approach 


components 




V,. =r CL«V a a* -\- CL~V a a~ 4- fl~\? * a~ 


V A = V A A 






+ Laao 1 #0 „ Iabo 


^AA_p = LAA_b +k 0 I_ AA0 + k.Q m LAB0 


where : #i = a 2 , a 0 — a , a n = 1 
a = exp(j2#73) 


wnere . fc 0 - # * 0m - 


Z — Z 

P 0 ==2M — =0m_ (for symmetrical lines: P 0 =l). 


c-g fault: 




^ _ imagiVju _ p [3a(/*A0 ~£oIabo)]*} 


(27c) 


/mag{(Z 1M / i4A _ p )[3a(/ AA0 -P 0 /aso)]*} 


where : 




Fault loop signals composed 


Fault loop signals composed 


in terms of symmetrical 


as in the classic approach 


components 








Z/lA_p = &iIaaI + ^2^AA2 + 






^AA_p = Z/U_c + I£qIaA0 ^^OiuLabO 


where : a x = a , a 2 =a 2 , a 0 =l 
a = exp(j2;r/3) 


where- it -^ola-^ila , Z 0 m 

—IIA 



P 0 = ~ r um (for symmetrical lines: P 0 = 1 ) . 

—OLA -^0m 
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10 



Distance to fault (d) for phase- to- phase- to ground faults can 
be calculated in two different ways, depending whether the 
pre- fault currents can be used or have to be avoided. 

1. Procedure for distance to fault calculation with the use of 
pre-fault measurements: 

a-b-g fault: 

imag{V^MilAAo -P 0 Labo)Y} 



imag{(Z l£A LAA_p)\KilAA0 -PqIabo)Y} 



(28a) 



where: W = =^ + ^ J 



AA2 



CALaa^DIaaz' 



A = l-a' 



B=l~a , C = l + a 2 , D = 1 + a 



Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as in the classic approach 


VAA^p ="iVaA1 +&2Y-AA2 +«0^AA0 


¥-AA_p -V-AA^a ~Y-AA_b 


, n —OLA r , ^Om r 

+ ry LAA0 1 ry —ABO 

±L\LA =-\LA 


!-AA_v = i-AA_a ~ I-AA_b j 


where: ^=1 — a 2 , a 2 = 1 — a , a 0 = 0 
a =exp(j2^-/3) 





to =~ M — = 22L (for symmetrical lines: £ 0 = 1) 

—OLA -Z-Om 



15 b-c-cr fault: 

imagiiZ^l^MlAAO -PoL ABO )Y} 



(28b) 
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A = a 2 -a , B = a-a 2 , C = -l, D = -l 



Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as in the classic approach 


¥-AA_p =£lKA41 +ZlY-AA2 






LAA_p - LAA_b ~LAA_c 


where: a x = a 2 -a t a 2 =a-a 2 , a 0 =0 
a =exp(j2;r/3) 





Z — Z 

£-0 =~ (for symmetrical lines: P 0 = 1 ) 



c-a-gr fault: 

^{(ZiiaLaa^MUaao -EoIabo)!} 
where : 



CAL^ + DLaai ' 
A = a-1 , B = a 2 -I, C = a + 1 , D = a 2 +1 



Fault loop signals composed 
in terms of symmetrical 
components 


Fault loop signals composed 
as in the classic approach 




¥-AA_p = ¥-AA_c ~~¥-AA_a 


Laa^p =£iZaai +& 2 LaA2 + 

, n —OLA j . 2om r 
— 1LA — 1LA 


-AA_p ~ I-AA_c ~I-AA__a 


where : a x = a - 1 , a 2 = a 2 - 1 , a 0 = 0 
a =exp(j2/r/3) 
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Eo = f 0LB f ° m (for symmetrical lines: £„ =1 > 



2. Procedure for distance to fault calculation without the use 
of pre-f ault measurements : 
a-h-g fault: 

d = imag[(V a +V b )(lAA0 -PqLab o)*} 



imag[Z iIA a a +L b +2k 0 I AAQ +2k 0m l AB0 )a AA0 -PqLabo)*] 



(29a) 



b-c-g fault: 

imag[Z l£A (l b +l c +2k 0 l MO +2k 0m lABo)aAAO ~P 0 Iabo)*1 



(29b) 



c-cL-g- fault: 

d = imag[(V c +V a XLAA0 -goiW*3 

imag[Z l£A a c +l a +2k 0 l AAO +2k Qm l AB o)a A Ao -RoIabo)*] 



where- Ic„ - — QLA ~ lLA Jr - — <>"* 
wrier e K 0 - , /c 0/?j - — 



It is also noted that while the above describes exemplifying 
embodiments of the invention, there are several variations and 
modifications which may be made to the disclosed solution 
without departing from the scope of the present invention as 
defined in the appended claims. 
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Tables 

TABLE 1. Coefficients for determining the fault loop voltage 
(V M p ) and current (Laa_p ) in terms of symmetrical components 
as defined in (2) and (3) . 



i Fault type 




« 2 




a-g 


1 


1 


1 j 


b-a 


2 

a 


CL 


1 


c-g 


a 


a 2 


1 


a-Jb, a-b-g 
a-b-c t a-b-c-g 


l-a 2 


l-a 


0 


b-c, b-c-g 


g 2 —a 


2 

a — a 


0 


c-a, c-a-g 


a-1 


a 2 -1 


0 


a = exp(j2#73) , j = V~T 



TABLE 1A. Fault loop voltage (VyU_ P ) current (Laa_ v ) 

expressed with using the classic approach. 



Fault type 






a-g 






| h-g 




I-AA_b + k-ol-AAO + komlABO 


c-g 




LaA_c +£o!/UO ^^OhiLabO 


a-b, a-b-g 
a-b-c, a-b- 
c-g 




LaA_o, —I-AA^b 


b-c, b-c-g 




I-AA_b ~ LaA__c 


c-a, c-a-g 




Laa_c ~Laa_o. 


AlLA ^ILA 

a=exp(j2^/3) , j = V r l 



TABLE 2 . Alternative sets of the weighting coefficients for 
determining a voltage drop across the fault path resistance 
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Fault 
type 


Set I 
( recommended ) 


Set: XX 


Set III 




Q.F2 


Q.F0 


&F1 


O.F2 


O.F0 


3.FI 


&F2 




a-er 


0 


3 


0 


3 


0 


0 


1,5 


1* 5 


0 


b-gr 


0 


3a '■ 


0 


3« 2 


0 


0 


1,5a 2 


1,5a 


0 


c-gr 


0 


3a 2 


0 


3a 


0 


0 


1,5a 


1,5a 2 


0 


II a ~ b 


0 


l-a 


0 


l-a 2 


0 


0 


0,5(l -a 2 ) 


0,5(1 -a) 


0 


! b ~° 


0 


2 

a— a 


0 


a 2 — a 


0 


o 


Ojig 2 -a) 


0,5(a -a 2 ) 


0 


c-a 


0 


a 2 -I 


0 


a-1 


0 


0 


0,5(a-l) 


0,5(a 2 -l) 


0 


I a-h-g 


l-a 2 


l-a 


0 


l-a 2 


l-a 


0 


l-a 2 


l-a 


0 


b-c-g 


2 

a —a 


2 

a— a 


0 


a 2 —a 


2 

a — a 


0 


2 

a —a 


a — a 2 


0 


c-a-g 


a-1 


a 2 — I 


0 


a -I 


a 2 - I 


0 


a-1 


a 2 -l 


0 


a-b-c-g 
(a-b-c) 


l-a 2 


0 


0 


l-a 2 


0 


0 


l-a 2 


0 


0 


a = exp(j2x 1 3) , j = V— T 



TABLE 3. Coefficients for determining a fault current 
distribution factor (6) 



SINGLE LINE (Fig.4) 


Z ?£oc 
— IAB 


£l = ~?L\L?L\AB ~~ (Z.\sA + )Zl\L 


Zliab "~ >oc 


£l = "~^lZf ^1 = 2lL + ^ IsB 

Mi =z lsA +z lsB +z lL 



PARALLEL LINES (Fig.5) 








Ll\ — Zl\uSZ-\sA +2l\ s B + £\LB&Ab) + Z-ILB&AbZ-IsBI 




Ml -Z.XLaZaLB&.AB +Il\Ia(Z±\sA + ZlUb) + £\LB&Ab(ZisA +Zl\sb) 




wnere. Ailb&ab — _ „ 



42 



WO 2004/001431 PCT/SE2003/001091 
Table 4 The recommended SET of the coefficients b Fl ,b F2 in 

relation to (2 6) 



T3*=an 1 f- 






a -sr 


0 


1 _ 


b-g 


0 


a 2 


c-g 


0 


a 


a - exp(j2^/3) , 



Table 5 The recommended SET of the coefficients a F1 , a F2 , a F0 
5 in relation to (26) 



I FAULT 


S.FI 


2F2 




a-gr 


0 


3 


0 


b-g 


0 


3a 


0 


c-g 


0 




0 


a — exp(j2;z73) , j = V-T 
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APPENDICES 

Al. DERIVATION OF THE COEFFICIENTS a } , a 2 , a 0 (TABLE 1) 
Single phase-to-ground fault: a-g fault 
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l-AA_v ~LAA_a + KqLaAO + koml-ABO = i-AA\ + LaA2 + IaAO + ~ 0lA _ ~~~ lLA LaAO + 'Z^ IaB0 = 

= LaM +LAA2 +=^lAAO +§ M -L AB 0=«dAAl + ^LaA2 +"0=^1^0 + *(>,„ f^ABO 
±l\LA —ILA —ILA 

Thus : a x = a 2 = a Q = 1 
5 Inter-phase faults: a-Jb, a-Jb-gr, a-Jb-c, a-jb-c-gr faults 

= (l - a 2 + (l - ajV^ = a 2 y ^ + a 2 + a 0 V M0 

^A4_p ~LAA_a -l~AA_b ~(LaA\ + IaA2 + IaAQ ) " ^ZaAI + £l>U2 +^AAo) = 
= & " ^ )^AA1 + (1 " ^UaA2 = filial + «2l/U2 + £o ^ I /WO 

Thus: a x = 1 - a 2 , a 2 =l-a, a 0 = 0 



10 A2 . Derivation of coefficients a Fl , a F2 , a F0 (TABLE 2) 

TABLE 2 contains three alternative sets (Set I , Set II, Set 
III) of the weighting coefficients, which are used for 
determining a voltage drop across a fault path. The 
coefficients are calculated from the boundary conditions - 

15 relevant for a particular fault type. This is distinctive that 
in all the sets the zero sequence is omitted (a F0 =0). This is 
advantage since the zero sequence impedance of a line is 
considered as the uncertain parameter. By setting a F0 =0 we 
limit adverse influence of the uncertainty with respect to the 

20 zero sequence impedance data upon the fault location accuracy. 
To be precise one has to note that this limitation is of 
course partial since it is related only to determining the 
voltage drop across a fault path. In contrast, while 
determining the voltage drop across a faulted line segment the 

25 zero sequence impedance of the line is used. 

a.-g fault, Figure 12: 

30 
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Taking into account that in the healthy phases: l F b =Lf_c =0 

what gives: l Fl =\(lF_a +<*I-F_b +a 2 I F _ c )=^{lF.a +£0 + 2 2 o) = -i/ F _ a 

Ln =\k F _ a + <k 2 L F _ b +<4-F_c)=\kp-a +« 2 0+ao)=|/ F _ a 

The sequence components are related: 1_ F} — L F2 = Lfo and 
finally: 

l F =L F _a = 3 L F 2 * thus: a Fl = 0 , a F2 =3 f a FQ =0 (as in the Set I 

from Table 2) 
or 

l F = I Fa =3/ F1 , thus: a F1 = 3, a F2 =0, a F0 =0 (as in the Set II 

from Table 2) 
or 

/ F =J F a =l,5/ F1 +l,5J F2/ thus: a Fl = 1,5 , a F2 =l,5, a F0 =0 (as in the 
Set III from Table 2) 

a-jb fault Figure 13a, 13b: 

The fault current can be expressed as: L F -L F ^a or as : 

Taking into account that in the healthy phase: L F c - Q and 
for the faulted phases: Lf _b = "If _a < what gives: 




The relation between Lfi and Lf2 is thus: 
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Finally: 

3 

thus: a F1 =0 , a F2 = l-a , a F0 = 0 (as in the Set I from Table 2) 



or 



5 l F =L F _^a = ^ a ^ Fi = (l - £ 2 Yfi 

thus: a Fl = l-a 2 , a F2 = 0 , a F0 =0 (as in the SET 2 from Table 2) 
or 

If =0,5/ F _ a +0,5/ F _ a =^i^/ F2 = 0,5(l-a)/ F2 + 0,5(l-a 2 )r F i 

thus: a F1 =0,5(l-a 2 ), a F2 = 0,5(l-a) , a F0 =0 (as in the Set III 
10 from Table 2) 

(a-b-g) fault, Figure 14: 



Lf -If_o ~i-F_b ~ Ufi + Lf2 + Lfo)-(o. 2 Ifi +0Lf2 + Ifo) = 

=(i-« 2 kFi+(i-«)ii 



-F2 

15 Thus: a Fl = l~a 2 , a F2 = l-a , a F0 =0 (as in the Sets I, II, III 
from Table 2) 

(a-b-c) or (a-b-c-g) symmetrical faults, Figure 15a, 15b, 15c : 

20 Taking the first two phases (a, Jb) for composing the voltage 
drop across a fault path one obtains : 

Lf = ^F_a ~I-F_b ~ (If\ +If2 + IfQ )~ ^Lf\ +QJ-F2 + Lfq)~ 

= (l-« 2 kFl+(l-^)lF2 
Thus : 

a F1 =l — a , a F2 = I — a , a FQ =0 

25 Additionally, if a fault is ideally symmetrical the positive 
sequence is the only component, which is present in the 
signals. Therefore, we have: 
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a Fl = 1-a 2 , a F2 = 0 , a FQ =0 (as in the Sets I, II, III from Table 
2) . 

A3 . Derivation of the complex coefficients in the fault 

5 current distribution factors for the positive (negative) 
sequence (Table 3) 

a) Case of the single line with extra link 45 between the 
substations (Figure 4) 

10 Let us determine the fault current distribution factor for the 
positive sequence (the fault current distribution factor for 
the negative sequence is the same) . The equivalent circuit of 
Figure 4 with indicated flow of currents for the incremental 
positive sequence is presented in Figure 16. 

15 

Considering the closed mesh containing: the local segment of 
the faulted line, the remote segment of the faulted line and 
the extra link between the substations, one can write down: 

dZ XL AI A1 +a-d)Z iL (AI A1 ~l F1 )-Z lAB Al AB1 = 0 

20 From the above equation the unknown current from the extra 
link between the substation can be determined as: 

Considering the closed mesh containing the source impedances 
(Z UA , Z UB ) and the extra link (Z lAB ) one can write down: 

25 Z lsA (AL M +AL ABi ) + Z lAB AL ABl +Z lsB (AI A1 +Al AB1 -7 F1 ) =0 

Substituting the previously determined unknown current from 
the extra link into the above equation one obtains : 
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~ F1 Lfi Mi 



where, as in TABLE 3 (Single line, Z lAB ) , we have: 

5 L\ = Z lL {z XsA + Z x JjB ) + Z lAB (z lL + Z ljB ) 
Mi =(z 1M +Z lj5 )(Z 1Afl +Z 1L ) + Z 1L Z 1AB 

If there is no extra link between the substations (Z 1AB — >oc) one 
has to consider the closed mesh containing the source 
impedances {Z lsA , Z lsB ) and both the segments of the faulted 
10 line [dZ 1L and (1 — d )Z lL ] . For this mesh one can write: 

(ZlsA +dZ 1L )AI A1 +\z lsB +(l-d)Z lL ]{AL A1 -I Fl ) = 0 

After rearrangements one obtains: 

k = AI -ai = K i d+L l 
~ F1 Lfi Mi 

where, as in TABLE 3 (Single line, Z 1AB — ) , we have: 
15 K X =-Z 1L 

Ki=z lsA +z lsA +z 1L 

b) Case of the parallel lines with extra link between the 
20 substation (Figure 5) 

Let us determine the fault current distribution factor for the 
positive sequence (the fault current distribution factor for 
the negative sequence is the same) . The equivalent circuit of 
25 parallel lines from Fig. 5 with indicated flow of currents for 
the incremental positive sequence is presented in Figure 17 . 
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The healthy parallel line (LB) and the extra link 55 (AB) , 

which are in parallel connection, have been substituted by the 

equivalent branch with the equivalent impedance: 



Z_\r.nZ_\ 

—1LB&AB - " . ry 

—ILB + ^1AB 



7 - ±L\LB*±IAB 

Zjy - - • " — 



5 Considering the closed mesh denoted by (AA, F, BA, BB, AB) one 
can write: 

^Z^ALaax + (1 - d)Z 1LA {ALaax -IfO- Z ILB8cAB AI WSlABI =0 

From the above equation the unknown current from the 
equivalent branch can be determined as : 

ZL1LB&AB 

Considering the closed mesh containing the source impedances 
(Z lsA , Z lsB ) and the equivalent branch (Z lLBScAB ) one can write 
down: 

^ua(^Laax + a Ilb&abO + Zilb&ab a Llb&abx +Z.isb( a Laai + a Ilb&abi ~Lf0 = 0 

15 Substituting the previously determined unknown current from 
the healthy line into the above equation one obtains: 

h _ 4Lm _ K l d + L l 

Lfi Mi 

where, as in TABLE 3 (Parallel lines), we have: 

^1 = ~£-ILa(Z-1sA +Z.UB +ZiLB&Ab) ' 

20 L x -Z lLA {Z lsA +Z ljB + Z lLB&AB ) + Z 1LB&AB Z lsBl 

Ki =Zi£aZilb&ab +Z lLA (Z lsA + Z lsB ) + Z 1LBScAB {Z lsA + Z lsB ) 
where: Z 1LB&AB = ^xlbZiab . 

In case if the extra link between the substations (Z lAB ) is 
not present one has to substitute : Z XLB&AB =Z lLB . 

25 
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